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I , Abstract. Calculations are presented which use the molecular R-matrix with pseudo- 

^ ' states (MRMPS) method to treat electron impact electron detachment and electronic 

O ■ excitation of the carbon dimer anion. Resonances are found above the ionisation 

I threshold of C2 with ^Hg and '^Ilg symmetry. These are shape resonances 

C/3 ' trapped by the effect of an attractive polarisation potential competing with a repulsive 

. ^ , Coulomb interaction. The Hg resonances are found to give structure in the detachment 

cross section similar to that observed experimentally. Both excitation and detachment 

(— I ' cross sections are found to be dominated by large impact parameter collisions whose 

O I contribution is modelled using the Born approximation. 

T— I ■ 

> : 

^ ■ 1. Introduction 

(N I 

^ J i The carbon dimer anion is unusual among small anions in having bound electronically 
^ i excited states. In addition electron scattering experiments performed in storage rings 
00 ! showed that was one of a number of diatomic anions which, at least temporarily, 
^ I are able to bind an additional electron leading to pronounced resonance structures in 
I their measured cross sections (Andersen et al. 1996, Pedersen et al. 1998, Pedersen 
^ I et al. 1999, Andersen et al. 2001, Collins et al. 2005). Similar results have also been 
5^ I found for systems containing the NO^ anion (?, ?). Given that the electron collisions 
occur against the background of a strongly repulsive Coulomb interaction, the occurrence 
of such resonance structures, which generally lie above the threshold for electron impact 
detachment, was unanticipated. 

Theoretically the treatment of electron collisions with these diatomic anions is 
complicated by the need to treat the region immediately above the electron detachment 
threshold. Such intermediate energy calculations are difficult because of the presence of 
two continuum electrons. Previous theoretical studies have used bound state methods 
either unadapted (Pedersen et al. 1999) or with an absorbing potential (Sommerfeld 
et al. 1997, Sommerfeld et al. 2000) to study the continuum states of Cl~ . 

In a preliminary publication (Halmova & Tennyson 2008) we reported calculations 
using the molecular R-matrix with Pseudo-states (MRMPS) method (Gorfinkiel & 
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Tennyson 2004, Gorfinkiel & Tennyson 2005) which found a number of low-lying 
resonances. In this work wc report fully on these calculations, inchiding details of 
the models we tested to produce stable calculations and to demonstrate that our results 
are robust. In addition we report results for electron impact excitation of 

2. Method 

The i?-matrix method is based on dividing coordinate space into two regions using a 
spherical boundary of radius a centred on the centre- of-mass of the target molecule 
(Burke & Bcrrington 1993, Burke & Tennyson 2005). The radius of the boundary is 
chosen so that the inner region contains all the electronic cloud of the target molecular 
states included in the calculation. 

Inside the i?-matrix sphere it is necessary to consider all short-range interactions 
between the N target electrons and the scattering one, such as exchange and electron 
correlation. In the outer region these effects are negligible so the scattering electron can 
be described as moving in the long-range multipole potential of the target molecule. 

The accuracy of this method is strongly dependent on the representation of the 
problem in the inner region (Tennyson 1996 &). In standard, low-energy calculations, 
the wave function for an 7V-I-1 electrons system in the inner region is given by the 
expansion: 

= ^^ajjfe$i(xi...xjv)wij(xjv+i) + ^bikXi{^i-^N+i) , (1) 

ij i 

where k represents the A;*'^ solution of the inner region Hamiltonian, A is the 
antisymmetrisation operator, Xj are the spatial and spin coordinates of electron i, 
Vjij are continuum orbitals (COs) which represent the scattering electron (Tennyson 
& Morgan 1999), aijk and hi^ are variational coefficients, <l>i is the wave function of the 
i^^ target state and Xi ^-re functions constructed from the target occupied and virtual 
molecular orbitals. These functions represent electron correlation and polarisation 
effects. In the first sum, the configuration state functions are constrained to give the 
correct (target) space and spin symmetry for the first A'"-electrons as well as the correct 
total, + 1 electron space-spin symmetry. Doing this requires special consideration of 
phase effects due to electron ordering in the wave function (Tennyson 1997). 

In the standard formulation of the R-matrix method all the solutions of the inner 
region problem are required to construct the R-matrix on the boundary (Burke & 
Berrington 1993). This presents a serious computational barrier for large calculations as 
diagonalising the entire Hamiltonian may not be feasible. For some of the calculations 
discussed below we used the partitioned R-matrix method (Tennyson 2004), which only 
requires the low-lying eigenvalues and which has been demonstrated to give good results 
for the electron - problem (Tennyson & Halmova 2007). 

Here we use the UK polyatomic i?-matrix code (Morgan et al. 1997) rather than the 
specialised. Slater orbital-based diatomic code. This is because the MRMPS method 
described below relies on the use Gaussian Type Orbitals (GTOs) and is therefore only 
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implemented in this code. The highest symmetry available in the polyatomic code is 
which is a subgroup of the true D^oh symmetry of C^. All calculations presented 
here were performed in symmetry. In the polyatomic suite target, continuum and 
MRMPS pseudo-continuum orbitals are all represented by a linear combination of GTOs. 
The target wave functions are expanded as a linear combination of the configurations 

$i(xi...xjv) = ^ Qfc0fc(xi...XAr) , (2) 

k 

where the Cik coefficients are determined by diagonalising the Hamiltonian of the 
molecular target. The quality of the target wave functions is dependent on the size of this 
expansion, as is indirectly, the quality of the scattering calculation (Tennyson 19966). 

The central idea of the MRMPS method is the augmentation of the close-coupling 
expansion (1) with extra "target" wave functions that represent pseudo-states. Unlike 
the usual target wave functions, these pseudo-states are not approximations to true 
eigenstates of the target, but are used to represent a discretised version of the electronic 
continuum. These pseudo-states arc obtained by diagonalising the target electronic 
Hamiltonian expressed in an appropriate basis of configurations, see below. 

As we are considering an anionic target, there are a finite number of bound target 
electronic states: three in the case Cj. This means that, unhke previous MRMPS 
studies (Gorfinkiel & Tennyson 2004, Gorfinkiel & Tennyson 2005), the pseudo-states 
are expected to all lie above the threshold to ionisation (electron detachment) where they 
represent the discretised target continuum. Of course, as the pseudo-continuum orbitals 
are added to the target basis they also influence the representation of the target states. 
It is a standard and tested assumption of the R-matrix with pseudo-states method 
that electron impact ionisation cross sections can be obtained by summing the cross 
sections for excitation of pseudo-states which lie above the vertical excitation threshold 
(Bartschat & Bray 1996). 

To generate configurations that describe an ionised target, the MRMPS method 
uses and extra set of orbitals called pseudo-continuum orbitals (PGOs). These orbitals 
are used to describe the ionised electron. The PGOs are expanded in terms an even- 
tempered basis set (Schmidt & Ruedenberg 1979) of GTOs centred at the centre of mass 
of the system. In this type of basis sets, the exponents of the GTOs follow: 

ai^aoP^'-'\ (3>1, i = l,....,L, (3) 

where by choosing different values of the parameters cuq and /3 different basis sets can 
be systematically generated. This is useful for checking the convergence and stability 
of the calculation and, of particular importance here, identifying physical resonances as 
distinct from the pseudo-resonances which are a known artifact of the RMPS procedure 
(Bartschat et al. 1996). 

Given that the GOs arc also represented by GTOs, to avoid problems with linear 
dependence it is necessary to remove those GO basis functions from the set whose 
exponents are greater than cuq. Even after this condition care has to be taken to ensure 
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that the three basis sets (target, CO and PCO) give a hnearly independent set of 
orbitals. To do this the PCOs are first Schmidt orthogonahsed to the target molecular 
orbitals (MOs) and then symmetric orthogonahsed among themselves. At this stage 
those orbitals with eigenvalues of the overlap matrix less than a deletion threshold 5, 
here taken as 4 x 10^^, were assumed to be linearly dependent and were removed from 
the basis. This procedure is repeated for the COs, here using 5 ~2 x 10~^. 

3. Calculations 

3.1. Target representation 

The starting point for the present calculations was our previous study of the electron - 
C2 system (Halmova et al. 2006). In that work the states of C2 were represented using 
the double-zeta plus polarisation (DZP) Gaussian basis set of Dunning (1970), natural 
orbitals and a complete active space configuration interaction (CAS-CI) in which the 
four Is electrons were frozen in the lag and lcr„ orbitals, and the remaining eight 
electrons were freely distributed among the 2o"c,, 3o"c,, 2cr„, 3o"„, Itt^ and Ivr^ orbitals 
giving configurations which can be written {lag lauY (^Cg 3(jg 2a"„ 3cr„ l7r„ Itt^)^. To 
treat C2 we added a PCO basis comprising 10 s, 10 p and 6 d orbitals. In this subsection 
all results presented used PCO exponents generated using q;o=0.17 and P—l.A. The 
calculations presented here are for C2 in its equilibrium geometry for which R — 2.396 oq 
(Huber & Herzberg 1979). 

The previous MRMPS studies considered electron impact ionisation of H2 and 
(Gorfinkiel & Tennyson 2004, Gorfinkiel & Tennyson 2005), which are both two electron 
systems. In this case the construction of ionised target plus PCO configurations is 
straightforward. This is not so here where there are many possible configurations that 
could be selected. The need to (a) get a good representation of the (pseudo)-continuum, 
(b) get good energies and wave functions for the physical states of the target, (c) obtain 
a balanced description between the target and scattering calculations and (d) keep the 
whole calculation computationally tractable meant that considerable experimentation 
was required. Table 1 summarises the models, giving the target configurations used in 
(2), which were tested as part of the present study. 

Besides choosing a basis set and a set of configurations, it is also necessary to build 
a set of target molecular orbitals. For standard scattering calculations these orbitals are 
normally ones associated with the A^-electron target under consideration. However for 
calculations involving ionisation the final target state only has N—1 electrons; experience 
(Gorfinkiel & Tennyson 2004), echoed by tests performed as part of this work, has shown 
that use of orbitals associated with the ionised target gives best results. We therefore 
tested two sets of C2 MOs, those generated from a self consistent field (SCF) calculation 
and natural orbitals (NOs) generated using the prescription given previously (Halmova 
et al. 2006). 

Table 2 gives energies for states of Cg calculated with the various models and 
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Table 1. Configurations used in the various target models tested. N is the size of the 
resulting Hamiltonian matrix for ^Ag symmetry and PCO means pseudo continuum 
orbital. 



Model 


N 


Configurations 


1 


1110 


(1^3 


lcr„)^ (2crg ScTg 4(Tg 2au 3au Ivr^ Itt^)^ 








2a, la„ 2a„ IttJ^^ (PCOs)i 


2 


14(1 




2(7, Idu 2au)^ ItVu^ {Sag 3au li^gY 








2ag lau 2(7^,)^ InJ {Sag 3a„ iTTg)^ 








2ag lau 2auf InJ (PCOs)i 








2ag lau 2(tJ8 {Sag Sa^ iTr,)^ (PCOs)^ 


(J 






lau 2(Tg)^ ItTu"^ (2(Tu Sag 3(j„ iTrg)^ 






(1^3 


lcr„ 2<Tp)^ IttJ (2cr„ Scr^ 3<T„ Itt^)^ 






(1^9 


lau 2agY InJ (2a„ Sa^ 3f7„ IttJ^ (PCOs)i 








la„ 2aj6 inj {2au Sag 3a„ IvrJ^ (PCOs)i 


4 


425 




lcr„)'^ (2cTg 3crg 2(7^, Sau Itt^ Itt^)^ 








2(7, la„ 2auf IttJ (PCOs)i 








2a, la„ 2au)' In J {Sag Sau Ing)' (PCOs)i 




597 




lauY {2ag Sag 2au Sau l^u l^gf 








lau) (2^3 2(j„ Itt^)'' (PCOs)^ 








lau)^ {2ag 2au IttJ^ (3(7, 3a„ IvrJ^ (PCOs)i 


6 


3111 


(1^9 


1(T„)^ (2(7, 3(7, 2(7u 3(7„ l7r„ Itt,)^ 








1(7„)4 (2(7, 2(7„ l7r„)« (PCOs)i 








1(T„)4 (2(7, 2(7„ (3(7, Sau l^gY (PCOs)i 








1(7„)4 (2(7, 2au (3(7, 3(7„ l7r,)2 (PCOs)i 


7 


20454 




lc7„ 2ag)^ (2a„ 3(7, 3cr„ Itt, Itt^)^ 








1(7„ 2(7,)6 {2au Sag Sau Itt, l7r„)6 (PCOs)^ 


8 


97500 




lau)^ (2(7, 3(7g 4(7, 2au Sau l^u l^g)^ 






(1C^9 


lauT (2t^9 Sa^ 4(7, 2(7„ 3(t„ Itt^ iTr,)^ (PCOs)i 



orbitals sets. These can be compared with the results of our previous study (Halmova 
et al. 2006) which, once nuclear motion effects had been taken into account, gave good 
agreement with the experimental results reported in Huber & Herzberg (1979). 

Considering each model in turn. Model 1 gave reasonable target energies but the 
pseudo-state configurations are very limited as they are obtain from a single electron 
in a PCO and a frozen target C2; this treatment is not consistent with the CAS-CI 
representation of Cg . Model 2 is more consistent in that 8 electrons are frozen in all 
configurations. However this much more limited model gave poor energies for the target 
states; indeed with SCF MOs it did not even predict the correct ground state for C^. 
Model 3 is built on Model 2 but with only 6 target electrons completely frozen; it gives 
rather a large scattering Hamiltonian and when tested did not give particularly good 
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Table 2. C2 ground state energies (negative numbers in E?i) and excitation energies 





(in ('\') for 


(lillort^iit mod 


("Is and ()rl)ita]s. 
















u 


J_> 
















1 


-75.61166 


- 


0.729 


- 


2.673 


- 


2 


1.951 


-75.53709 


-75.56659 


0.850 


7.955 


5.813 


6 


-75.59480 


-75.58782 


A CT AT 

0.597 


0.611 


3.115 


2.549 


4 


-75.60956 


-75.60944 


0.468 


0.419 


2.734 


2.506 


5 


-75.62675 




0.883 




2.717 




6 


-75.65777 




0.912 




2.908 




7 


0.123 


-75.71419 


-75.66055 


0.686 


1.823 


2.621 


8 


-75.72075 




0.689 




2.621 




a 




-75.67213 




0.557 




2.355 



" Previous calculations (Halmova et al. 2006). 



eigenphase sums. Model 4 was built as a hybrid between Models 1 and 2 using the 
target configurations from the former and pseudo-state configurations from the latter. 
This model gave good target energies and a relatively small Hamiltonian size for the 
scattering problem. This model formed the basis of our preliminary study (Halmova & 
Tennyson 2008) and became our workhorse for test calculations. Model 5 is a slightly 
enlarged version of model 4 and behaves similarly. 

Use of the partitioned R-matrix method allowed us to explore target models which 
implied significantly larger Hamiltonians for the scattering problem. We therefore tested 
the effect of gradually expanding the CAS-CI used to generate the pseudo-states. The 
limit of the process is Model 8, in which the same extended CAS is used in both parts 
of the calculation. However, its corresponding scattering model is far too big to be 
tractable (see table 4) and was not pursued. Our attempts to construct intermediate 
models, gives numbers 6 and 7. Model 7 did not give good target energies. We therefore 
decided to concentrate on use of models 4, 5 and 6. 

As can be seen from table 2, the use of C2 SCF MOs gives results of similar 
quality to NOs, in contrast to more usual calculations were NOs give significantly better 
results. This is perhaps not surprising since the NOs were constructed to give a good 
representation of a range of electronically excited states of C2 which is not our purpose 
here. The calculations reported below use SCF MOs unless otherwise stated. 

There is one further target property which proved to be of considerable importance 
for the scattering calculations, that is the long range polarisability of C2 target. 
Gorfinkicl & Tennyson (2004) showed that the MRMPS method converges the long- 
range target polarisability in a way that standard close-coupling expansions do not. As 
is discussed below, a good representation of the polarisability is essential to get a good 
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Table 3. Isotropic polarisabilities of C2 


in ao for different models, orbitals and PCO 




basis determined by {ao,(3). 




Model 


ao=0.17 (3=1A q;o=0.15 P=1A 


q;o=0.17 /3=1.5 q;o=0.17 /3=1.3 




SCF MOs C2 NOs SCF MOs 


r-H -| — ^ Ti /%- y^ /— y y^T — 1 jr y^ "TV -T- y^ 

SCF MOs SCF MOs C2 NOs 


1 


26.42 




2 


28.58 109.51 




3 


67.49 12.24 




4 


32.24 32.48 


31.24 31.98 18.80 


5 


25.00 




6 


24.55 





physical model of the electron - molecular anion collision as it provides the dominant 
attractive term in the interaction. 

Table 3 shows the isotropic polarisabilities predicted by various of our target models. 
We could find no literature value for this parameter. Indeed it is not straightforward to 
calculate it with standard electronic structure codes, we tried, as they are not generally 
set up to treat molecular anions. The table shows that the calculated polarisability 
is generally stable with changes to the PCO basis but not with other aspects of the 
target model. In part this is because both excited states of C2 are dipole allowed from 
the ground state and low-lying. This makes the calculated polarisability particularly 
sensitive to the transition dipole and precise excitation energies of these states. Our 
results suggest that the true polarisabilty of C^ at i? = 2.396 ao is between 25 and 32 
ag, of this less than 10 ag arises from coupling to the two physical electronically excited 
states of C^. 

So far we have concentrated on our calculations of the physical properties of the 
target. There is however one other property of importance for the scattering runs: the 
spectral coverage of the continuum by the pseudo-states. Figure 1 compares energy 
levels generated by various versions of Model 4. The behaviour shown here, that the 
energy levels are stable to choice of PCO basis but very sensitive to the target model 
used, was also shown in the other comparisons we made (Halmova 2008). A common 
feature of all the comparisons is the sparsity of states directly above ionisation. In 
principle one could get pseudo-states in this energy region by significantly expanding 
the R-matrix box size and the associated PCO basis; such a calculation was not deemed 
computationally tractable at present However, as can be seen from the results presented 
below, our pseudo-state distributions lead to very small near-threshold electron impact 
detachment cross sections, which appears to be in agreement with the observations. 
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a„=0.n |3=1.4 



a„=().15 (3=1.4 



a=0.n p=1.5 



a=0.n p=1.3 



a|=().17 |3=1.3 
NOs 




B E 



A"n 



IP 



Figure 1. Target state distribution for C2 for Model 4 for various (ao,/3) values as 
indicated in the figure. Energies are relative to the X ground state of C2 . SCF 
MOs of C2 were used except where otherwise indicated. The horizontal line marked 
IP denotes the ionisation threshold. 



3.2. Scattering model 

Test calculations using Model 1 for R-matrix spheres of radius a = 10 ao and 13 ao 
showed that a = 10 ao gave stable results and this was used for all further calculations. 
CO's were taken from Faure et al. (2002) with the largest exponent functions removed. 
This basis contains functions with £ up to 4 ie g orbitals. Use of two different sets of 
PCOs were tested: 

4-140^ 2-762. 2-763« l-46ip 4-86i„ 2-5hg 2-5b2g la„ 
4-240^ 2-1062. 2-IO63. l-7hg 4-126i„ 2-863^ 2-862^ la„ 

in D2h notation. The eigenphase calculated with the larger PCO set was slightly higher, 
so this was used for further calculation. 

To go from the target models detailed above to the inner region R-matrix wave 
functions of (1) we followed the prescription developed previously (Tennyson 19966) for 
standard R-matrix calculations. This prescription, which has been used successfully 
in many studies, is designed to provide a balanced treatment between the target, A^- 
electron, and scattering, + 1-electron problems. Table 4 details the configurations 
generated for the 1-electron calculations associated with each of the models described 
above. Also given is the size of the final Hamiltonian, a crucial parameter in determining 
the tractability of the calculation. Although our Hamiltonian construction algorithm 
is very efficient meaning the computational demands of this step of the calculation 
only depend weakly on the number of target states (and pseudo-states) included in 
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the expansion (Tennyson 1996a), the diagonahsation step still represents a major 
bottleneck. As a result, the largest previously published molecular R-matrix calculation 
was restricted to TV ~ 28000 (Rozum et al. 2003). 

Model 8 gives an extremely large final step Hamiltonians and no attempt was 
made to actually solve for it. Many calculations were performed for Models 1 to 6 for 
testing purposes. Calculations using the partitioned R-matrix method (Tennyson 2004) 
exphcitly considered the lowest 1000 solutions which were found to be sufficient to span 
scattering energies up to about 40 eV. 

Scattering calculations were performed for both singlet and triplet symmetries. 
After some experimentation it was decided to include 114 target states in the final 
close-coupling expansion. This number is evaluated in symmetry so counts states 
which are degenerate in Dooh symmetry twice. 66 of these states are doublets of which 
the lowest 3 are the physical states of C^. The 48 quartet states only couple to triplet 
calculations. These states span energies up to about 19 eV above the target ground 
state for both Models 4, 5 and 6. 

Outer region calculations were performed by propagating the R-matrices to 100 ao 
and then matching to asymptotic Coulomb functions. This large number of states 
makes the outer calculations slow compared to standard R-matrix calculations but we 
were still able to generate results for sufficient energies to map out the various resonance 
structures, both real and pseudo. 

3.3. Born Approximation 

The interaction between the scattering electron and the C^ anion is long-range and 
repulsive. It can therefore be expected that many of the collisions will occur with a 
large impact parameter. Such collisions are not well represented in the calculations 
described above, all of which use a CO basis set and hence partial wave expansion 
truncated at £ = 4. However collisions with high values of electron angular momentum 
(£ > 4) are amenable to a much simpler treatment since in these collisions the scattering 
electron can be assumed not to penetrate the charge cloud of the C2 target. Under these 
circumstances the collision cross sections are determined by purely long-range effects. 

To allow the important effect of collisions with £ > 4 in the inelastic cross sections 
presented here, we used a simple top-up procedure based on the Born approximation 
(Baluja et al. 2001). There are a number of ways applying such a Born top-up procedure 
but since C2 is non-polar, meaning rotational effects will be small, and the non-resonant 
cross sections at low I are very small, meaning that interference effects for these ^'s will 
not be important, the simplest procedure based on the direct augmentation of the cross 
sections was used. This procedure was applied to electron impact electronic excitation 
of all physical and pseudo states which are dipole allowed, that is to all excitations to 
states of and ^Ily symmetry. Only transition dipole moments were considered for 
the long-range potential. For the pseudo-states this extra excitation cross sections were 
taken as supplementing the electron impact detachment cross sections. This increase 
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Table 4. Configurations used for the various scattering models tested. N is the 
size of the resulting Hamiltonian matrix with ^Ag symmetry and PCO means pseudo 
continuum orbital and CO means continuum orbital. 



Model 



N Configurations 



21705 



6447 



52270 



6575 



12283 



12283 



823823 



too big 



lag IfTu)'* (2o-g Sag iag 2au 3au Ittu IvTg)^ (COs)-*- 

lag 2ag lau 2cr„ l7r„)i2 (PCOs)l (COs)i 
lag IcTu)^ (2o-g Sag iag 2au 3o"u ItTu ItTj)!" 

lag 2ag lau 2(t„ ItTuY"^ (PCOs)^ 

lag 2ag lau 'iau l7ru)12 {-^^^ 4^^^ 3^^ (pcOs)l 

lag 2a g lau 2auf Itt^^ (3o-g 3a„ iTTg)! (COs)l 

lo-g 20-9 lau 2auf liTu-^ {Sag Sa^ liTgf (COs)l 

lo-g 20-9 lau 2auf liTu* (PCOs)l (COs)l 

lag 2ag lau 2cr„)S Itt^S (30-5 Sau li^gY (PCOs)l (COs)l 

lag 2ag lau 2cr„)S ItTu* {Sag Sctu IJTg)^ 

lag 2ag lau 2(t„)« Itt^* [Sag Sau iTTgY (PCOs)i 

ICa 2ag lau 2fT„* l7r„* (PCOs)^ 

2ag lau 2<T„)8 Ircu^ {Sag Sau l-JVgf (PCOs)i 

2ag lau 2(T„)» l7r„3 (3^^ Sau iTTg)! (PCOs)^ 



la, 
la, 

lag lau 2(Tg)® ItTu* {2cr„ 



^ ^ l7rg)3 (COs)i 

lO-g lO-u 20-g)6 iTVu^ (2<T„ 3<Tg 3(T„ iTTg)* (COs) ^ 

la, ' 

IcTg lau 2ag)^ ItTu"' (2(T„ So-g Sau iTTg)* 



1<T„ 2(Tg)6 l7r„4 {2au Sag Sau iTVgf (PCOs)l (COs)l 
lau 2agf iTTu^ {2au Sag Sau iTTgf (PCOs)i (COs)l 



lag lau 2ag)^ ItTu* (2(7„ 3o-g Sa, 



IWgf {PC0S)1 



lau 2agf' l7r„4 (2a„ Sag Sau liTg)'^ (PCOs)2 
lau 2agf l7r„^ (2(7„ Sag Sau iTVgf (PCOs)l 
lau 2agf l7r„3 (2ct„ Sag Sau iTTgf {PCOs)^ 
Ictu)** (2(7g Sag 2au Sau ItTu l^rg)^ (COs)^ 
2<T9 lau 2auf Ittu'^ (PCOs)l {COs)l 
2ag lau 2auf lvr„^ {Sag Sau IVgY (PCOs)l (COs)l 
Ictu)** (2o-g Sag 2au Sau li^u iTg)"'" 
2ag lau 2auf ItTu'^ (PCOs)^ 
2<Tg lau 2auf ItTu^^ {Sag Sau iTVgf (PCOs)l 
2<T3 lau 2auf 1-kJ {Sag Sau iT^gY (PCOs)^ 
1(T„)4 {2ag Sag 2au Sau li^u l-^g)'^ (COs)i 
1(T„)4 (2CTg 2au I'Kuf (PCOs)i (COs)i 
1(T„)4 (2^9 2au liTuf (3<Tg Sau It^g)^ (PCOs)i 
lauY (2a'g Sag 2cr„ Sau l7r„ iTTg)^" 
lauf {2ag 2au iT^uf (PCOs)^ 
lau)^ {2ag 2au liVuf {Sag Sau iTrg)^ (PCOs)l 
lau)^ {2ag 2au ItTuY {Sag Sau liTg)^ (PCOs)^ 
Ictu)** (2o-g 3o-g 2au Sau lvr„ iTrg)'^ (COs)"' 
lau)^ {2ag 2au Itt^)* {PCOs)1 (COs)l 



(COs)i 



ICTg 
ICTg 
ICTg 
ICTg 
lO-g 
lO-g 
lo-g 
ICTg 
lo-g 
lo-g 
lo-g 

lag 
lag 

lag 
lag 
lag 
lag 
lag 

lag „, „ ~, s - - , ^ - - , 

lag lau)"" {2ag 2(7„ liTuf {Sag Sau li^gY (PCOs)l (COs)l 
lag lauf {2ag 2au liVuf {Sag Sau iTTgf (PCOs)l (COs)l 
lo-g lau)"^ {2a g Sag 2au Sau ItTu l'^g)^° 
lag lauf {2ag 2(t„ l7r„)« (PCOs)^ 
l<Tg lau)* {"iag 2au liruV {Sag Sau li^gf (PCOs)i 
iCTg lau)* (2(Tg 2au iTTu)"^ {Sag Sau liTg)^ (PCOs)^ 
lo-g lo„ 2og)S (2a-„ 3o-g Sctu iTTg Ittu)'^ (COs)l 
lOg lo„ 2og)6 {2au Sag Sau iT^g It^u)^ (PCOs)i (COs)i 
lo-g lou 2og)^ {2au Sag Sau l^g Ittu)* 
lo-g lo„ 2og)6 (2o-„ 3o-g 3o„ iTTg Ittu)^ (PCOs)i 
lOg lo„ 2og)6 {2au 3ag Sau iTtg iTTuY' (PCOs)2 

lOg lo„)* (20-9 30-9 20u 3o„ ItTu iTTg)^'' 

lOg lau)* {2ag 2au l7r„)» (PCOs)^ 

lOg lau)* (2og 2o„ iTTu)'^ {Sag Sau Ivrg)^ (PCOs)l 
lOg lau)* (2og 2o„ iTTu)'^ {Sag Sau lvrg)i (PCOs)2 
lo-g lou)** (2o-g 3o-g 4CTg 2au Sau Ittu iTTg)^ (COs)l 
log lau)* (2og 3og 4og 2o„ 3o„ iTTu l7rg)S (PCOs)l (COs)l 

lOg lOu)'* (20g 30g 40g 20u 30u ItT^ iTTg)-"' 

log lo„)'' (20-9 30-9 4og 2o„ 3o„ liTu iTTg)" (PCOs)i 
lo-g lau)* {2ag Sag 4og 2o„ 3o„ l7r„ iTrg)^ (PCOs)^ 
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Table 5. Energy, E,., and width, T, both in eV, of resonance states of : our results 

are l)ase(l calculatious witli n()=0.17 and /i=1.4. 



Previous work Model 4 Model 5 Model 6 



Syimuotrv" 


Er 


r 




Syumiotry 


Er 


r 


Er 


r 


Er 


r 




3.5 


0.3 


a 




4.86 


0.65 


4.75 


0.57 


4.91 


0.67 




10.0 


2.1 


b 




10.92 


0.52 


10.59 


0.57 


11.46 


0.82 




10.0 


3-4 


b 




9.71 


1.14 


9.54 


0.99 


10.44 


1.43 



Theory (Sommerfcld ct al. 2000). 
^ Experiment (Andersen et al. 1996, Pedersen et al. 1998, Pedersen et al. 1999). 



turns out to be substantial showing that the majority of this process does in fact occur 
through long-range collisions. 

4. Results 

A very large number of different calculations were attempted. For space reasons only 
only those performed using Models 4, 5 and 6 will be reported in detail. As a major aim 
of this study was to characterise the resonances observed in the storage ring experiments, 
it is this aspect of the work we discuss first. We will then give our results for total cross 
sections for both electron impact electronic excitation and electron detachment. 

4.1. Resonances 

All models tested showed a resonance feature at about 5 eV. However only the 
more sophisticated models showed stable resonances in the energy region probed by the 
storage ring experiments. 

Figure 2 gives a series of Model 4 calculations for ^83^ symmetry contribution to 
the total electron impact electron detachment cross section as a function of energy. Both 
singlet and triplet contributions are dominated by a series of resonance features: several 
narrow resonances which show a strong dependence on the PCO basis set parameters 
and a single broad resonance in each case. In contrast to the narrow resonances, the 
position and width of the broad resonance features is stable with respect to the choice 
of PCO basis; indeed even the calculation with q;o=0.17 and /3=1.5 which shows a much 
enhanced resonance structure, a known signature of a poorly converged PCO basis 
(Bartschat et al. 1996), gives a similar resonance to the other calculations. Figure 3 
gives a similar comparison for the total electron impact detachment cross section where 
again both the singlet and triplet resonance features are clearly visible. The poorer 
quahty of the calculation based on the use of C2 NOs can clearly be seen. 

Inspection of plots of eigenphase sums (not given), showed one further clear 
resonance at energies just above the ionisation threshold of symmetry. Parameters 
for this resonance and the ones near 10 eV were obtained by fitting their eigenphases to 



Electron collisions with 



12 



T \ ^ \ I \ ^ \ ^ \ ^ \ ^ ^ 




Energy (eV) 

Figure 2. Partial electron impact detachment cross sections of C2 for Model 4 for 
final states of ^63^ (^Hg) total symmetry with different values of {ao,P); ionisation to 
singlet and triplet final states are shown to illustrate their distinct resonance structures. 
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u 



0.5 



ag=0A7 p=1.4 
a„=0.15 P=1.4 
a„=0.17p=1.5 
a„=0.17p=1.3 
a^=0.17p=1.3 CjNOs 



Energy (eV) 



Figure 3. Total electron impact detachment cross sections of for Model 4 with 
different values of (ao,/3)- 



a Breit-Wigner formula (Tennyson & Noble 1984). The resulting resonance parameters 
are given for Model 4 in table 5, where they are compared to previous work. 

Sommerfeld et al. (2000) performed careful absorbing potential calculations which 
found a single resonance close the ionisation threshold of • We predict a similar 
resonance, albeit at slightly higher energy; this agreement can be considered satisfactory 
given the lower level of configuration interaction included in our study. However we 
find two further resonances not found by Sommerfeld et al. (2000). These resonances 
lie close in energy to those detected experimentally (Andersen et al. 1996, Pedersen 
et al. 1998, Pedersen et al. 1999). In fact the experiments give somewhat different 
resonance characteristics depending on whether they measure the detachment cross 
section, which gives a width of 2.1 eV, or the smaller dissociation cross section, which 
gives widths between 3 and 4 eV. However to really test whether these resonances agree 
it is necessary to make a more direct comparison with the experiments which means 
comparing with the measured cross sections. 



4-2. Cross sections 

Figure 4 gives a comparison between our calculations, performed with three different 
models, and the measurements. What is immediately apparent is that without the Born 
correction our cross sections are very significantly lower than the measured ones. This 
supports the assertion made above that the majority of the ionisation occurs through 
large impact parameter collisions. 
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Energy (eV) 

Figure 4. C2 electron impact ionisation cross sections. Experimental data due to 
Pedersen et al. (1999) is given by circles (without dissociative channels) and squares 
(with dissociative channels); the dashed lines represent our cross sections without Born 
correction; the solid line represent our cross sections with Born correction. The arrow 
indicates the location of the ionisation threshold. Top panel: Model 4 with q:o=0.17 
and (3=1A; Middle panel: Model 5 with q:o=0.17 and /3=1.4; Bottom panel: Model 6 
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All the models shown give reasonable agreement with the measurements of Pedersen 
et al. (1999). For completeness we give the measured results with and without the 
consideration of simultaneous ionisation and dissociation. Our calculations implicitly 
give the cross section for all processes that involve ionisation and therefore should 
account for both channels. As the dissociation channel is relatively unfavoured, these two 
measurements are close together, the difference is certainly smaller than the accuracy 
of our calculations. 

It can be seen that all three models give total ionisation cross sections in good 
agreement with the measurements. In this context it should be noted that the cross 
sections should all go to zero below the ionisation threshold denoted by the arrow. 
Furthermore our resonance structures are very similar to, if a little higher in energy 
than, the observed ones. For reasons discussed for the lower-lying resonance, we 
would expect our resonance energies to be a little too high. This aspect of the agreement 
can therefore also be considered to be very good. 

Finally figure 5 gives total electron impact electronic excitation cross sections for 
the two bound states of C2 . These excitations are both dipole allowed and the cross 
sections are again completely dominated by the high impact parameter, dipole-driven 
transitions which we characterise with the Born approximation. Indeed as we predict 
essentially no resonance structure in these cross sections, a simple Born calculation 
should suffice for this problem. 

5. Discussion and conclusions 

Our calculations detect the clear signature of 3 resonances: one of symmetry just 
above the ionisation threshold and resonances with singlet and triplet spin and 11^ 
spatial symmetry at about 10 eV. Since they all lie above the excited electron states 
of C2 , these must be regarded as shape resonances. The standard model for shape 
resonances in electron - molecule collisions relies on a centrifugal barrier to provide 
the trapping potential. Such resonances have two characteristics: they are not found 
for s-wave scattering and they occur, often at too high an energy, in simple static 
exchange calculations. The present resonances are absent from the test static exchange 
calculations we ran and the resonance involves s-wave scattering. 

All this implies a modified picture of a shape resonance. The situation here is that 
the scattering occurs under the influence of a strongly repulsive Coulomb potential. The 
well required to (temporarily) trap the scattering electron should therefore be the result 
of some local attraction. Our calculations suggest that this is indeed what happens 
and the attractive component is provided by polarisation interactions. This is why 
characterising the resonances require such a sophisticated representation of the target 
wave functions. 

Although representation of the resonances requires a very sophisticated treatment 
of the target, this is not the situation for the cross sections. Both the electron impact 
electronic excitation and electron detachment cross sections are dominated by long range 
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Figure 5. Electron impact electronic excitation cross sections for excitation to the 
A ^n„ and B states without (upper) and with (lower) Born correction; note the 
very different cross section scales for the two panels. Calculations are for Model 4 with 
ao=0.17 and (3=1.3. 
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collisions which can be modelled using a simple Born approximation. It should however 
be noted that our treatment of the electron impact electron detachment process does 
involve the use of a discretised continuum via the MRMPS method which is then used 
to provide the necessary transition dipoles for the Born calculation. 

In summary our calculations on electron impact electron detachment of 
reproduce both the resonance structure and cross sections observed over a number of 
years in storage rings (Andersen et al. 1996, Pedersen et al. 1998, Pedersen et al. 1999) for 
this process. These calculations are the first to achieve either of these goals. To do this 
we had to combine a sophisticated treatment of the C2 target including representation 
of target continuum states using the molecular R-matrix with pseudo states (MRMPS) 
method with a method for calculating the effects of high impact parameter collisions. 
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